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A  COST  OPERATOR  APPROACH  TO  MULTISTAGE 
LOCATION- ALLOCAT ION 
by 

Robert  V.  Nagelhout  and  Gerald  L.  Thompson 


The  multistage  factory-warehouse  locat lon-al locat Ion  problem  Is  to 
decide  on  locations  of  warehouses  and  shipping  amounts  from  the  factories 
through  the  warehouses  to  meet  customer  demands  in  such  a  way  that  the 
total  fixed  plus  variable  costs  are  minimized.  Capacity  constraints  at 
factories  and  warehouses  are  also  Imposed. 

We  present  a  cost  operator  algorithm  for  solving  this  problem.  The 
algorithm  takes  into  account  the  network  structure  and  the  submodularity  of 
the  objective  function.  Computat ional  results  with  problems  taken  from  the 


literature  as  well  as  new  problems  are  provided. 

* 


1  .  lnt  roduc  t  Ion 


In  this  paper  we  pre»cnt  a  cost  operator  algorithm  for  aolvlng  multistage 
factory  warehouse  locat lon-al location  problems.  The  decision  variables  cor* 
respond  to  the  warehouse  locations  and  the  shipping  amounts  from  the  factories 
through  the  warehouses  and  Into  the  demand  centers.  The  problem  Is  to  minimize 
the  total  fixed  costs  of  locating  warehouses  plus  the  totsl  variable  shipping 
cost,  subject  to  possible  capacity  restrictions  at  the  factories  and  the  ware¬ 
houses,  plus  the  demand  requirements  at  the  demand  centers. 

The  algorithm  takes  advantage  of  the  network  structure  of  the  supply 
and  demand  constraints  and  the  submodularity  of  the  objective  function.  We 
use  cost  operators  [22).  to  facilitate  the  movement  up  and  down  the  search 
tree.  This  greatly  reduces  the  amount  of  time  spent  solving  transportation 
subproblems,  which  frequently  comprises  up  to  90T  of  the  computational  burden. 

In  Section  2  we  give  a  problem  formulation  and  we  discuss  some  of  the 
more  recent  research  related  to  the  multistage  location  problem.  In  Section  3 
we  point  out  lower  bounds  and  fathoming  rules  obtainable  from  submodular  set 
functions.  In  Section  4  we  transform  the  multistage  location  problem  Into  a 
transportation  problem  format,  and  we  show  how  cost  operators  can  be  used  to 
generate  feasible  solutions.  Section  5  contains  a  description  of  the  cost 
operator  branch  and  bound  algorithm.  In  Section  6  we  give  an  example,  and 
in  Section  7  we  provide  extensive  computat Iona  1  experience  on  probleats  from 
the  literature  plus  some  multistage  problems  of  our  own. 

2 .  Problem  formulation 

We  describe  the  location  problem  to  be  studied  In  this  paper  as  well  as 
similar  models  which  have  been  presented  In  the  literature.  We  use  the  following 


notation : 
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q  •  number  of  factories;  I*  ■  ll,...,q) 

t  •  number  of  warehouses;  I  •  { q  + 1 . m’ 


r  "  number  of  demand  centers 

m  ■  q  ♦  i  *  1 ;  n  *  t  +  r  +  1 ;  J '  ■  : J  •  {<*•!,.. 


fixed  cost  of  opening  warehouse  i 


cost  per  distance 
per  unit  of 
shipping  from 


,  factory  1  to 
'warehouse  J 

factory  1  to 
'  demand  center  J 

warehouse  1  to 
demand  center  J 


id' 

j«j' 

id' 

j«J 

id 

jcJ 


.n} 


Aj  •  capacity  of  factory  1 

S.  •  capacity  of  warehouse  i 

d (  -  demand  at  location  J 

•  amount  shipped  from  location  i  to  location  J 

1  if  warehouse  1  is  open 

y  • 

1  0  otherwise 

The  multistage  or  intermediate  location  problem,  which  we  will  call 
Problem  P.  can  be  formulated  as: 


Minimize  Z(T)  •  Z  Z  c  x  ♦  I 

Til  Id'LT  Jej'-J  J  let 

Vi 

(1) 

subject  to 

id' 

(2) 

jL  •5‘”‘  i0 

icT 

(3) 

jd 
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-  ,  *  1  x14  *  0  t«T  (5) 

he  I  hl  jej  1J 

xtJ  >  0  id'  U  I,  jej'  U  J  (6) 

1  UT  <7\ 

^1  •>  0  otherwise 

Constraints  (2)  and  (3)  ensure  that  the  amount  shipped  out  of  a  factory  or  a 
warehouse  should  not  exceed  its  capacity.  Constraints  (4)  require  that  the 
demand  at  each  demand  center  be  satisfied.  Constraints  (5)  are  the  standard 
“conservation  of  flow"  constraints  which  require  that  the  amount  shipped  Into 
each  warehouse  equals  the  amount  shipped  out.  Constraints  (6)  and  (7)  are  the 
nonnegativity  and  Integrality  constraints  respectively.  The  objective  In 
problem  P  Is  to  minimize  the  total  transportation  costs  from  shipping  plus 
the  total  fixed  costs  of  opening  warehouses,  while  satisfying  the  customer 
demand.  When 


S,  <  I  d  (8) 

J«J  J 

we  say  that  warehouse  1  is  capac ltated  since  It  cannot  satisfy  all  of  the 
demand  by  itself.  The  factories  can  also  be  either  capacitated  or  uncapacitated 
depending  upon  the  size  of  A ^ .  Note  that  In  the  multistage  formulation  when 
I  *  9,  the  total  warehouse  capacity  need  not  exceed  the  total  demand  since 
units  can  also  be  shipped  directly  from  the  factories  to  the  demand  centers. 

A  considerable  amount  of  research  has  been  performed  on  uncapacitated, 
capacitated,  and  mixed  (partly  capacitated)  location  problems.  It  seems  that 
each  problem  which  Is  studied,  depending  upon  the  objective  function  and  the 
capacity  assumption,  exhibits  its  own  charscterlstlcs  and  yields  a  different 
algorithmic  approach.  Thus  the  literature  on  the  many  different  location 
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problems  la  enormous.  In  this  section  we  will  discuss  recent  literature  only 
for  those  problems  which  are  closely  related  to  P. 

Research  on  multistage  location  problems  has  been  limited.  Geoffrion 
and  Graves  (12]  described  and  tested  an  algorithm  for  solving  multistage 
multicommodity  distribution  systems  using  Benders'  decomposition.  Ellwln 
and  Cray  (8)  described  and  tested  an  algorithm  for  solving  single  stage 
d'  »  0)  location  problems  and  proposed,  but  did  not  test,  an  algorithm  for 
solving  multistage  problems.  Host  of  the  papers  In  the  literature  have  con* 
centrated  on  the  single  stage  version  of  P  for  which  I*  •  0.  This  has 
been  called  the  simple  (capacitated  or  ut.^apac ltated )  warehouse  location 
problem  which  we  denote  by  PI.  Some  of  the  models  which  have  been  studied 
impose  added  configuration  constraints  on  PI  which  restrict  the  total 
number,  and  different  combinations  of  warehouses  which  can  be  opened.  For 
example,  a  constraint  frequently  used  is: 

I  y.  <  K  (9) 

id 

where  K  Is  an  upper  limit  on  the  total  number  of  warehouses  which  can  be 
opened.  In  the  special  case  where  •  0  for  Id  and  each  warehouse  Is 
uncapacitated,  PI  with  constraint  anmended  Is  called  the  K  medn 

problem.  Both  heuristic  [14]  and  exact  methods  (5],  (16],  (21],  have  been 
proposed  for  solving  the  K  median  problem. 

Earlier  attempts  at  solving  the  uncapacitated  PI  concentrated  upon 
the  relaxation,  Pl ' ,  of  PI  where  (7)  is  replaced  by: 

0  <  yt  <  l  Id  (10) 

(see  [6],  (7|).  The  basic  Idea  was  to  solve  PI  by  Imbedding  Pl'  Into  an 


1  »«im— 
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lmplielt  enumeration  scheme.  Since  then  Pi'  has  been  called  the  weak  linear 
programming  relaxation  of  PI.  It  is  termed  "weak"  because  there  is  another 
linear  programming  relaxation  of  PI  which  has  proven  to  be  much  stronger  than 
PI  .  To  describe  the  latter,  we  add  to  problem  PI  the  constraints: 


<  min  •s1*djJy1 


lc I,  JcJ 


(11) 


Let  us  denote  by  Pi  the  problem  in  which  the  constraints  (11)  are  added  to 
PI  .  For  uncapac it ated  location  problems  Pi'  is  called  the  strong  linear 
programming  relaxation  of  PI.  (Note  that  for  uncapacitated  location  problems 
alniS^.d^J  •  mln’*,li  ■  1).  Pi'  is  stronger  than  Pi'  in  the  sense  that  the 
gap  between  the  optimal  values  of  PI  and  Pi''  Is  normally  much  smaller  than 
the  gap  between  the  optimal  values  ot  PI  and  Pi'.  Also,  It  Is  often  the 
case  for  uncapacltated  location  problems,  that  a  solution  to  Pi"  will  satisfy, 
or  almost  satisfy  (7).  That  is,  after  solving  Pi"  almost  all  of  the  fixed 
charge  variables  will  be  naturally  Integer.  However,  even  though  optimal 
solutions  to  PI  tend  to  be  close  approximations  to  an  optimal  solution  for 
PI,  researchers  have  attempted  to  avoid  solving  Pi'  directly,  because  it 
has  an  ercrmous  number  of  constraints  (there  are  mxn  constraints  of  the  type 
(ll>).  Schrage  (20),  has  proposed  a  method  for  solving  linear  programs  which 
handles  constraints  of  the  type  (11)  implicitly,  thereby  reducing  storage 
requirements.  Still  the  time  required  to  solve  Pi"  by  the  simplex  method 
can  be  excessive  and  could  cause  difficulties  In  an  Implicit  enumeration  scheme.. 

Instead  of  solving  Pi"  directly,  heuristic  methods  have  been  proposed 

which  find  feasible  solutions  to  the  dual  of  Pi"  .  To  describe  these  approaches 

* 

let  us  denote  by  a  feasible,  and  by  X^j  an  optimal  solution  to  problem  PI. 
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L« c  Z(X)  represent  the  objective  function  value  for  solution  X.  Also  let 
0  be  the  dual  problem  to  problem  Pi"  .  Then  by  duality  theory  we  know  that: 

)  -  Z(xJ)  >  ZlXp)  (12) 

* 

where  X^  is  an  optimal  and  X^  a  feasible  solution  to  D.  Thus  a  feasible 
solution  to  D  provides  a  valid  lower  bound  on  the  optimal  objective  function 
value  of  PI.  Subgradient  (a)  and  dual  ascent  methods  [3],  (9)  have  been  used 
to  find  good  feasible  solutions  to  D.  In  the  case  where  PI  Is  uncapac ltated 
one  can  easily  compute  a  low  cost  primal  solution  to  PI  after  having  found  the 

feasible  solution  to  D.  Furthermore  It  Is  often  the  case  when  PI  is 

a  * 

uncapac ltated  that  the  gap,  Z(Xp^)  *  Z(X^),  Is  very  small  or  even  zero, 
which  makes  the  additional  work  to  get  optimal  solutions  small.  Thus  the 
relaxation  PI  has  been  very  effective  In  solving  the  simple  uncapacitated 
warehouse  location  problems. 

Several  algorithms  also  exist  for  solving  capacitated  location  problems. 

Aklnr  and  Khumawala  111.  proposed  and  tested  an  implicit  enumeration  algorithm 

which  uses  Pi'  as  a  relaxation.  Ellwln  and  Cray  ( 8]  described  and  tested 

a  branch  and  bound  algorithm  which  uses  duality  properties  of  PI  and  bounds 

obtainable  from  the  suboodular lty  of  the  objective  function  (1)  for  fathoming. 

Culgnard  and  Spielberg  (13]  have  generalized  the  dual  ascent  method  of  Bllde* 

Krarup  (3)  and  Erlenkotter  (9]  to  the  capacitated  version  of  PI.  They  use  a 

relaxation  very  similar  to  Pi"  which  contains  the  constraints  (11).  They 

solved  some  randomly  generated  problems  and  found  that  the  zero  gap  phenomenon 
*  * 

between  Z(Xpj)  and  Z(Xp)  occurs  less  frequently  with  capacitated  or  mixed 
problems,  than  It  does  with  uncapacitated  problems.  Other  relaxations  besides 
Pi'  and  Pi"  have  also  been  used  to  solve  PI.  Geoffrlon  and  McBride  [11] 
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have  considered  a  location  model  for  which  PI  Is  a  special  case  and  have 
used  a  Lagranglan  relaxation  combined  with  Implicit  enumeration  to  solve  It. 
Recently  Nauss  (17)  used  a  Lagranglan  relaxation  to  solve  PI  with  excellent 
computational  results. 

In  the  following  sections  we  will  describe  an  algorithm  for  solving  P. 

The  general  approach  Is  similar  to  the  one  used  In  (8)  to  solve  PI.  We  do 
not  solve  a  relaxation  of  P.  Instead  we  make  use  of  some  lower  bounds  obtain¬ 
able  from  the  submodular  property  of  2,  together  with  some  other  fathoming 
rules,  to  enumerate  explicitly  a  subset  of  the  solutions  to  P.  The  movement 
In  the  search  tree  from  one  solution  to  another  Is  facilitated  by  applying  cost 
operators  (21)  to  P.  This  significantly  reduces  the  amount  of  effort  required 
to  solve  P,  since  most  researchers  have  found  that  the  majority  of  the  time 
Involved  In  solving  capacitated  location  problems  Is  spent  solving  transportation 
subproblems.  Ellwln  and  Cray  [8]  have  tested  and  shown  that  In  many  cases  over 
9 OX  of  the  time  required  to  solve  a  test  problem  Is  spent  solving  transportation 
subproblems.  In  the  next  section  we  discuss  some  lower  bounds  on  the  value  of  Z, 
some  of  which  are  utilized  In  the  algorithm  described  In  Section  5. 

3 .  Oblectlve  Function  Lover  Bounds 

Many  of  the  lower  bound  properties  which  we  present  here  have  been  dis¬ 
cussed  In  (2)  and  (10),  in  the  context  of  the  simple  uncapacitated  location 
problem.  These  lower  bounds  and  their  properties  are  also  useful  In  solving 
other  kinds  of  mathematical  programming  problems. 

There  are  many  properties  which  can  be  used  to  define  submodular  set 
functions.  For  a  discussion  of  them  see  (18).  To  define  one  such  function 
let  Z  be  a  real  valued  function  defined  on  the  finite  set  of  subsets  of  I. 
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For  notatlonal  convenience  we  define: 

iZ.(A)  -  Z(AUB)  -  Z  (A)  AC  I,  B  -  1  (13) 

o 

Then  any  of  the  following  properties  can  be  used  to  define  a  submodular 
function: 

(l)  Z(AJB)  ♦  ZCA^B)  >  Z(A)  ♦  Z(B)  for  all  subset#  A.BCI. 

(11)  AZ^A)  <  AZ  (B>  VACBCI,  lcl  CIA ) 

(111)  Let  {b. . B^j  be  a  partition  of  B-A;  then 

r 

Z(A)  >  Z(B)  -  I  iZ.  (B  -  B.  )  A-B.l.  (15) 

k-1  ®k  k 


In  most  expositions  property  (1)  Is  taken  to  be  the  definition  of  submodularlty , 
and  properties  (11)  and  (111)  are  shown  to  be  equivalent.  Details  are  omitted. 


In  the  context  of  this  paper,  Z(A)  represents  (1).  The  fact  that  Z(A) 
Is  submodular  was  proved  In  [15).  Note  that  in  P,  Z(0)  represents  the  value 
of  the  solution  where  all  of  the  shipments  originate  from  the  factories.  When 
l’  •  0,  as  In  PI,  one  must  be  careful  in  defining  Z(A)  If  A  happens  to 


be  Infeasible  set  of  warehouses,  l.e.,  when 


I  S  < 

i€A 


dj.  In  this  case  we 


we  define  a  "dummy  factory"  which  is  alwavs  available  to  service  the  demand 


centers  but  at  a  high  shipping  cost.  This  makes  Z(A)  large  enough  so  that 
the  value  of  Z (A)  for  any  infeasible  A  Is  at  least  as  large  as  the  value  of 


any  feasible  solution.  The  dummy  factory  approach  preserves  the  submodularlty 
of  Z  and  also  yields  a  different  value  of  Z  for  solutions  having  different 
degrees  of  lnfeas ibi 1 ltv .  As  we  will  see  later,  this  is  helpful  In  deciding 


which  warehouses  to  open  when  we  are  working  with  Infeasible  sets  of  warehouses. 

Properties  (14)  and  (IS)  can  be  Interpreted  as  adding  or  subtracting  ware¬ 
houses  from  a  given  set  A  of  open  warehouses.  For  example  property  (14) 


says  Chat  the  addition  of  warehouse  i  to  the  set  A  decreases  the  total 
cost  by  at  least  as  much  as  the  addition  of  warehouse  1  to  the  set  B  when 
A  -  & .  Thus  it  is  similar  to  the  "decreasing  returns  to  scale"  condition  in 
economics . 

Properties  (lu)  and  (13)  can  be  used  to  characterize  solutions  to  P 
and  to  derive  lover  bounds  on  the  value  of  any  solution  to  P.  The  algorithm 
to  be  described  in  Section  3  searches  for  sets  Til  which  have  the  following 
two  properties: 

(1)  uZt<T)  2  0  Vic  1-T 

(11)  AZt<T  -  (i}>  <  0  V  i  c  T 

Sets  T  m  I  whivh  satisfy  (i)  and  (ii)  are  such  that  the  addition  to  T  or  the 
deletion  from  T  of  a  single  facility  does  not  cause  Z  to  decrease.  Clearly 
any  optimal  solution  to  P  satisfies  (1)  and  (11).  In  a  sense  properties  (1) 
and  (ii)  characterize  the  set  of  all  "locally  optimal"  solutions  to  P.  The 

globally  optimal  solution  is  the  best  locally  optimal  solution,  which  must  be 

found  by  a  search  process.  One  of  the  factors  which  make  it  hard  to  find  the 
globally  optimal  solution  is  that  there  are  many  locally  optimal  solutions  which 
have  nearly  optimal  objective  function  values.  Thus  in  any  enumeration  procedure, 
a  considerable  amount  of  effort  is  normally  required  to  eliminate  these  nearly 
optimal  solutions  from  consideration.  In  a  practical  sense  however,  it  may 
be  true  that  the  nearly  optimal  solutions,  say  those  within  IT  of  optimality 
may  Indeed  be  as  valuable  or  "as  optimal"  to  a  decision  maker  as  a  globally 
optimal  solution.  Given  the  inaccuracies  in  the  cost  data  and  other  environ¬ 
mental  and  political  factors  which  must  be  taken  into  consideration,  it  would  be 
desirable  to  have  not  only  an  optimal  solution  to  P,  but  in  addition  a  list  of 
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solutlons  which  are  nearly  optimal.  We  will  point  out  how  this  may  be  ac¬ 
complished  using  the  algorithm  of  Section  5. 

In  the  initialization  phase  of  any  location  algorithm,  two  rules  can 
be  applied  in  order  to  permanently  open  or  close  warehouses. 

Rl'LE  1  ■  If  uZ^(I  -  1 1  * )  3  for  any  1  «  I  then  the  warehouse 

1  will  be  open  in  some  optimal  solution  to  P. 

To  see  this  suppose  that  T  is  an  optimal  set  of  warehouses  and  i  t  T.  By 
property  (li) 

iZ^T)  -  Z(T  -  { i  j )  -  Z(T)  i  iZt(l  -  { 1 } )  £  0, 

thus  T  U  d  ;  is  also  an  optimal  solution. 

Rl'LE  2  ■  If  IZ^d  )  _•  0  for  any  1  C  I,  where  1  is  the  set 
of  warehouses  opened  by  application  of  RULE  1,  then  warehouse  1 
will  be  closed  in  some  optimal  solution  to  PI. 

The  Justification  for  RULE  2  is  similar  to  that  for  RULE  1.  We  will  show, 

In  the  next  section  how  the  testing  of  RULE  1  and  RULE  2  requires  only 
the  application  of  two  cost  operators  for  each  warehouse.  In  many  cases,  as 
will  be  seen  in  Section  7,  RULES  1  and  2  can  be  used  to  fix  open  and  closed  a 
large  portion  of  the  warehouses  in  an  optimal  solution. 

Another  property  of  suhmoduiar  set  functions  which  can  be  derived  from 

(li)  and  (15)  is  the  following:  Let  \Aj,....Arl  and  (Qj . )  be 

partitions  of  A  -  T  such  that  -  A (  for  each  1  •  l,...,t  and  some 
J  •  1 , . . .  ,r .  Then 

r  t 

Z(T)  >  Z (A )  -  I  £Z.  (A-A, )  >  Z (A)  -  Z  £Z.  (A-Q.)  (16) 

k-l  Ak  k  "  k-1  Ok  k 


for  all  subsets  T  satisfying  T  -  A  ~  I.  The  quantity  on  the  right  hand  side 
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of  (IS)  provide*  a  valid  lower  bound  on  the  value  of  Z(T),  for  any  T  »  I . 
Property  (6)  say*  when  we  u«e  a  more  refined  partition  of  A  -  T,  we  get  a 
weaker  lower  bound.  Suppose  we  let  A  •  I.  Then  the  most  refined  partition  of 

1  -  T  Is  clearly  i  { l  ^  } . Utl)  where  lfc  «  I  -  T  and  t  •  |l  -  T|  . 

Then  for  this  partition  (16)  yields: 

t 

Z  (T)  >  Z(I)  -  I  (I  -  UJ)  .  (17) 

k-1  1 i  k 

Among  the  class  of  lower  bounds  (15).  (17)  Is  the  weakest.  Notice  that  once 

Rule  1  has  been  applied  all  of  the  terms  on  the  right  hand  side  of  (17)  have 

been  calculated.  Thus  (17)  can  be  applied  at  any  time  atter  Rule  1  has  been 

tested  without  any  extra  computational  effort.  A  stronger  lower  bound  than  (17) 

could  be  obtained  with  some  added  computat lonal  effort  by  partitioning  I  Into 

sets,  A^,  of  sire  two  and  calculating  iZ^  (I  -  A^).  We  have  not  yet  tested 

k 

this  idea. 

Another  fathoming  device  which  provides  an  upper  bound  on  the  maximum 

number  of  warehouses  In  an  optimal  solution  to  P  can  be  obtained  as  follows: 

i*i  * 

Let  t »  1 1  -  I  | .  Let  f  ,...,f  ,  k  c  I  •  I  ,  be  a  nondecreasing  ordering 

1  *  * 

of  the  fixed  charges  not  In  I  .  Define 


DV 


-  AZ  (1 
J1  J1 


<Ji»  *  fl*  Ji 


«  I  -  I 


(18) 


and  let  DV  . DV  be  nonincreasing  ordering  of  the  DV  .  (Note  that 


J 


1 


Jt 


J. 


DV  <  0.)  The  quantity  -DV  ,  represents  the  smallest  possible  Incremental 

1  ”  u 

savings  on  the  total  shipping  cost  when  warehouse  Is  opened.  Let  Z  be 

any  upper  bound  on  the  optimal  objective  function  value  for  P.  Then  an  upper 

* 

bound  on  the  number  of  warehouses  In  an  optimal  solution  to  P  (containing  I  ) 

i  *i  * 

is  1 1  |  ♦  k 


where , 
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*  *  *  * 
k  t-k  k  +1  t-k  -1 

Z(I)-  I  t.  ♦  I  f.  ♦  E  lDV(  |<ZU<Z(l)  -  II,*  I  f.  ♦  I  |dv  ! 
it  l  '  l m  1  *1  s«l  *  t«I  t»l  *4  s«  1  Js 

Note  that  the  first  two  t«rai  on  the  right  hand  side  ol  (19)  give  the  smallest 

* 

possible  shipping  cost,  the  next  term  is  the  smallest  sum  of  k  •*  1  fixed 

charges,  and  the  last  term  is  the  smallest  possible  cost  of  closing  down 

* 

t-(k  ♦  l)  warehouses.  Therefore  the  right  hand  side  of  (19)  is  a  lower  bound 

i*t  * 

on  the  value  of  an  optimal  solution  to  P  containing  I  ♦  k  ♦  1  warehouses. 

u 

It  this  lower  bound  exceeds  Z  then  any  optimal  solution  will  contain  at  most 

*.  * 

I  |  k  open  warehouses. 

In  the  algorithm  of  Section  5  we  use  as  fathoming  devices  bounds  obtained 
from  (17)  on  the  value  of  an  optimal  solution,  and  bounds  obtained  from  (19)  on 
the  number  of  warehouses  in  an  optimal  solution. 

a.  ot  Cost  Operators  to  Solve  Problem  P 

Clven  the  choice  of  a  subset  T  -  1,  problem  P  becomes  an  ordinary 
transshipment  problem.  Rather  than  resolving  this  problem  each  time  T  changes, 
we  use  the  operator  theory  of  parametric  programming  122)  to  change  the  problem 
and  derive  the  new  optimal  solution  s lmu 1 t aneous ly ■ 

We  shall  say  that  cell  (1,))  has  been  tlxed  out  of  the  basis  when  a  cell 
cost  operator  has  been  applied  to  the  problem  and  its  solution  so  that  the 
cost  has  been  driven  to  where  M  Is  so  large  that  x^  •  0  In 

any  optimal  solution  to  the  new  problem. 

We  also  say  that  cell  (1,J)  has  been  fixed  In  the  basis  when  a  cell 
cost  operator  has  been  applied  to  the  problem  and  its  solution  so  that  the 
cost  Cj,  has  been  driven  to  -M,  where  M  Is  so  large  that  x^  •  Mln(S^.d^) 

In  any  optimal  solution  to  the  new  problem. 


(19 
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F Igure  l  shows  an  example  of  a  2',3«4  multistage  location  problem 
Rows  l  and  2  correspond  to  the  factories,  rows  3,  4,  and  3  to  the  warehouses, 
and  row  b  represents  the  dummy  factory.  Columns  1,  2  jnd  3  correspond  to  the 
warehouses,  columns  4  through  8  to  the  demand  centers,  and  column  9  is  a  slack 
column.  Rows  1  and  2  contain  the  costs  of  shipping  from  the  factories  to  the 
warehouses  and  from  the  factories  directly  to  the  demand  centers.  Cells  (3,1), 
(4,2),  and  (5,3)  contain  the  fixed  charge  of  the  corresponding  warehouse, 
divided  by  its  capacity.  In  any  solution  to  P,  these  cells  will  contain  the 
unused  warehouse  capacity.  Cells  (3,2),  (3,3),  (4,1),  (4,3),  (5,1)  and  (5,2) 
cannot  be  used  because  of  their  large  costs;  in  effect  these  arcs  have  been 
removed  from  the  problem.  The  cells  in  rows  3  through  5,  and  columns  4  through 
8,  contain  the  costs  of  shipping  from  each  of  the  warehouses  to  each  of  the 
demand  centers  plus  the  proportional  fixed  charges.  Notice  that  some  of  the 
cells  In  dummy  factory  row  b  and  slack  column  9  contain  two  costs.  This  can 
be  explained  in  the  following  manner.  To  obtain  the  solution  where: 


-  0 

t 

c  •  -M 

in 

and 

Cm,l-q 

-  1 

set 

c  -  M 

In 

and 

c  . 
m,  l-q 

(In  Figure  l,  m  *  b,  n  •  8,  q  »  2,  and  1  «  (3,4,5}). 

To  see  this  consider  the  problem  shown  in  Figure  2.  By  solving  the 

transportation  problem  In  Figure  2  we  would  obtain  an  optimal  solution  to  P 

when  T  «  (2,3),  i.«.,  when  warehouse  1  is  closed  and  warehouses  2  and  3  are 

open.  Notice  that  call  (3,8)  has  a  cost  of  -M,  thus  In  the  optimal  solution 

n  .  >  S  ,  which  has  the  effect  of  closing  down  warehouse  1.  Cell  (6.1)  has 
3"  3 

a  cost  of  -M  which  causes  all  of  the  demand  In  column  1  to  he  satisfied  by 
the  d ■  xnmv  factory  i.e.,  xftj  ■  S^.  Thus  setting  c^g  •  -M  and  c^j  ■  -M 


•In¬ 


effectively  closes  down  warehouse  1.  On  the  other  hand  for  warehouse  2, 
c ,  _  •  M  and  c  ,  •  M,  so  that  x  „  ■  x.,  •  0.  This  causes  the  demand  in 

column  2  to  be  satisfied  from  rows  1,  2,  or  4.  The  shipments  from  rows  1  and  2 

represent  units  being  shipped  from  tactories  1  and  2  to  warehouse  2.  The 
shipment  from  row  4  In  colvimn  2  represents  the  unused  warehouse  capacity,  and 

it  is  charged  at  the  proportional  fixed  charge  rate.  Notice  that  if  cell  (4,2) 

contains  x, ,  units  then  exactly  $.  -  x, ,  units  can  be  used  to  ship  from 

*Z  4  *4  a» 

warehouse  2  to  the  demand  centers.  This  is  exactly  the  amount  which  is  made 
available  to  warehouse  2  from  the  factories.  Thus  the  conservation  of  flow 
equations  (3)  have  been  satisfied.  Because  a  proportional  amount  of  the  fixed 
charge  is  assigned  to  both  the  used  and  unused  parts  of  the  warehouse  capacity, 
the  total  fixed  charge  Is  covered  in  any  feasible  solution. 

Figure  5  contains  an  optimal  solution  to  the  2*4*3  examples  solved  in 
Section  6.  In  Figure  5,  T  •  { 1 , 3 ;  so  that  warehouses  1  and  3  are  opened,  and 
warehouses  2,  4,  and  5  are  closed.  Factory  1  is  not  used  at  all  so  that  x^  ^  ■ 
Factory  2  ships  21  units  to  warehouse  1,  40  units  to  warehouse  3,  and  lb  units  dir 
ectly  to  demand  center  1.  Warehouse  1  ships  all  21  of  the  units  it  received  from 
Factory  2  to  demand  center  3.  The  remaining  unused  4  units  of  capacity  at  ware¬ 
house  l  are  in  cell  (3,1)  and  are  charged  at  the  proportional  fixed  charge  rate. 
All  of  the  flow  for  warehouses  2,  4,  and  3  is  in  the  last  column,  and  their  demand 
is  satisfied  entirely  by  the  dummy  factory.  Warehouse  3  is  used  to  capacity 
shipping  22  units  to  demand  center  2,  and  18  units  to  demand  center  4.  In  this 
example  each  of  the  demand  centers  is  supplied  from  a  single  warehouse.  This  is 
not  the  case  In  general. 

The  idea  behind  the  cost  operator  approach  to  solving  P(T)  Is  to  open 


the  warehouses  in  T  and  close  those  in  I  -  T  by  fixing  in  or  out  of  the 
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basis  each  of  the  cells  (m,l)  through  (m, 4)  and  (q  +  1.  n)  through  q ,  n). 

For  example  the  problem  P(I),  In  which  all  of  the  warehouses  are  open,  can  be 
obtained  by  solving  the  transportation  problem  shown  In  Figure  3.  Then  to 
apply  Rule  1  In  Section  3  we  use  two  cost  operators  for  each  1  <  1.  For 
example  we  would  calculate  iZ^(I  -  (3))  •  Z(I)  -  Z(I  -  (.  3 ; )  by  fixing  In 
cells  (3,8)  and  (6,1).  Doing  this  would  yield  a  solution  to  the  problem 
shown  In  Figure  2,  and  permits  the  evaluation  of  Rule  1.  In  general,  the 
addition  of  a  warehouse  to  a  given  set  or  the  deletion  of  a  warehouse  from  a 
given  set  requires  the  fixing  In  or  out  of  two  cells.  The  amount  of  work  needed 
to  fix  two  cells  In  or  out  using  cost  operators  Is  much  less  than  the  compu¬ 
tational  effort  of  solving  a  transshipment  problem  from  scratch.  Since  the 
branch  and  bound  search  algorithm  to  be  described  In  the  next  section  requires 
the  solution  of  transshipment  problems  for  many  sets  T,  the  total  computational 
effort  saved  by  the  fixing  In  and  fixing  out  procedure  Is  very  large. 

Also  we  should  mention  that  in  the  case  where  I*  •  9,  as  In  single 
stage  problems,  only  one  cost  operator  is  needed  to  open  or  close  a  warehouse. 

In  the  single  stage  formulation  columns  1  through  L  and  rows  1  through  2  are 
absent.  Thus  In  order  to  fix  In  or  out  warehouse  1  we  need  only  apply  a 
cost  operator  tc  cell  (l.n). 

5 .  The  Cost  Operator  Algorithm. 

To  describe  the  cost  operator  branch  and  bound  algorithm  we  use  the  fol¬ 
lowing  notation: 

l  ■  level  of  the  search  tree 

Q  *  set  of  open  warehouses 

List(4)  •  list  of  warehouses  which  may  be  opened  on  level  l 
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1 

u 

Z  •  current  upper  bound 

X  ■  current  beet  solution.  1 

Ue  begin  by  setting  4  •  1.  Q  -  t  and  Z°  -  ZC9)  which  is  the  (high) 

cost  of  supplying  ell  demends  from  the  fectorles  directly.  Then  Rule  1  (see 

* 

Section  3)  Is  epplled  end  we  let  Q  ■  I  .  Then  for  eech  1  «  l-Q  we  spply 
the  lover  bound  test  (17)  by  setting  T  •  Q  -  (.  1 } ,  end  checking  to  see  whether 

ZU  >  Z(l)  -  E  AZ.(I  -  (k]).  (20) 

kcI-T 

If  the  right  hend  side  of  (20)  (which  Is  e  lower  bound  on  the  value  of 

ft  u 

Z (Q  w  U )))  exceeds  Z  ,  then  warehouse  1  can  be  permanently  closed.  For 

each  1  for  which  (20)  holds  true  we  apply  cost  operators  to  calculate  iZ^(Q). 

If  2iZj(Q)  _  0  then  we  permanently  close  warehouse  1.  Then  we  let  List(l)  be 

the  set  of  all  1  such  that  AZ^Q)  <  1.  (We  place  the  warehouses  In  List  (l) 

In  order  of  non- ir.c reas lng  objective  function  values.)  Next  we  remove  the 

last  warehouse.  1,  from  List  (l)  and  replace  Q  by  Q  U  {l}.  We  let  ZU  ■  Z (Q ) , 

update  X-  If  List  (1)  Is  empty  then  we  stop.  Otherwise  we  calculate  k*  as  In 
* 

(19).  If  k  •  l  then  we  stop.  Otherwise  we  replace  4  by  l  ♦  1  and 
continue  to  the  next  level.  At  each  level  after  the  first,  we  perform  the  lower 
bound  test  (20)  for  each  warehouse,  1.  In  List  (4.-1)  by  letting  T  ■  Q  'J  (.1 } . 

If  warehouse  1  passes  the  test  (20)  we  calculate  ^Z^(Q).  Again  we  let 

List  (4)  be  the  set  of  warehouses  such  that  2>Z^(Q)  <0  In  nonincrcaslng  order. 

If  List  (4)  Is  not  empty  we  remove  1,  the  last  element  of  List  (4);  let 
ZU  •  Z(Q  'J  (i } )  and  update  X.  If  Llat  (4)  la  now  empty  we  backtrack  by  re¬ 
placing  4  by  4-1,  removing  the  last  warehouse  paced  In  Q  and  then  continuing 

*  * 

as  before.  Otherwise  we  calculate  k  as  In  (19).  If  4  >  k  we  bracktrack 
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••  Just  described.  Otherwise  we  replace  t  by  4+1 ;  replace  Q  by  Q  U  ( i } , 
and  continue  as  before.  The  algorithm  terminates  when  List  (1)  has  been  emptied. 

Now  we  formally  state  the  cost  operator  multi-stage  locatlon-al location 
algorithm 

Step  (1).  (Initialisation).  Set  4  •  1 ;  rU  •  Z(0).  Apply  Rule  1.  Let  Q  •  I  . 

For  each  1  c  I-Q  apply  test  (20),  letting  T  -  Q  U  {l}.  For  each  1 

which  satisfies  (20)  calculate  iZ^lQ).  Order  those  1  such  that 
iZ^(Q)  •'  0  In  nonlncreas Ing  order  and  place  them  In  List  (1).  If 
List  (1)  empty  go  to  Step  (5).  Otherwise  let  k  be  the  last  warehouse 
in  List  (1).  Let  ZU  •  Z(Q  L‘  {k});  update  X.  Remove  warehouse  k  from 
List  (1).  If  List  (1)  Is  empty  go  to  Step  (5).  Otherwise  go  to  Step  (2). 

*  *  * 

Step  (2).  (Bound  by  k  ).  Calculate  k  as  In  (19).  If  4  >  k  go  to  Step  (6). 

Otherwise  replace  4  by  4-fl;  replace  Q  by  Q  w  (k)  and  go  to 

Step  (3). 

Step  (3).  (Forward  branching).  For  each  warehouse  1  In  List  (4-1)  apply 
test  (20)  letting  T  ■  Q  J  (i).  For  each  1  which  satisfies  (20) 
calculate  iZ^( Q).  List  those  1  such  that  Z.Z^(Q)  <0  In  non-in¬ 
creasing  order  and  place  them  In  List  (4).  If  List  (4)  is  empty  go 
to  Step  (6).  Otherwise  let  k  be  the  last  warehouse  In  List  (4). 

Let  z'*  "  Z(Q  -  Ik});  update  X.  Remove  warehouse  k  from  List  (4). 

If  List  (4)  Is  empty  go  to  Step  (4).  Otherwise  go  to  Step  (2). 

Step  (4).  (Backtracking).  Replace  4  by  4-1.  Remove  from  Q  the  last 

warehouse  placed  In  Q.  If  4  ■  1  and  List  (1)  contains  only  one  ware¬ 
house  go  to  Step  (5).  Otherwise  let  k  be  the  last  warehouse  In  List  (4). 
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Remove  k  from  List  (<■).  If  List  (>-)  empty  go  to  Step  (4).  Otherwise 
replace  Q  by  Q  U  {k};  replace  (.  by  4+l  and  go  to  Step  (3). 

Step  (3).  (Termination)  Stop.  The  current  solution  X,  Is  optimal. 


Notice  that  each  time  Step  (3)  Is  performed  a  new  list  of  feasible 

solutions  to  P  Is  available.  Many  times  these  solutions  have  nearly  optimal 

objective  function  values  and  thus  they  might  be  worthwhile  to  save. 

* 

After  calculating  I  ,  the  cost  operator  algorithm  proceeds  to  open,  one 
at  a  time,  those  warehouses  which  cause  the  largest  decrease  in  the  objective 
function  value.  This  Is  continued  until  at  level  K ,  no  further  decrease  In 
the  objective  function  value  Is  possible  (l.e.  List  ((.)  Is  empty).  Then  we  move 
to  the  backtracking  Step  (4).  The  best  feasible  solution  obtained  by  the  cost 
operator  algorithm  before  the  first  backtracking  is  called  the  greedy  solution. 
We  denote  by  Z&  the  objective  function  value  of  the  greedy  solution.  In  (18), 
a  worst  case  bound  on  z’  was  derived  for  submodular  set  functions  which  Is, 


* 


Z  <0 )-  Z 


♦  k  -  1 


t  *i  * 

ll  I  ♦  k 


(21) 


♦  * 
where  Z  is  the  optimal  objective  function  value  and  k  Is  obtained  from  (19). 

In  practice  however,  the  actual  percentage  error  obtained  by  the  greedy  solution  is 

much  smaller  than  the  worst  case  bound.  In  Section  7  we  will  see  that  In  most 

cases  the  greedy  solution  value  Is  within  .3  percent  of  optimality. 

Finally  we  should  point  out  that  it  would  be  trivial  to  "reverse"  the 
cost  operator  algorithm  so  that  Instead  of  starting  with  all  of  the  warehouses 
closed,  and  opening  them  one  at  a  time,  we  could  start  with  all  of  the  warehouses 
open  and  close  them  one  at  a  time.  This  could  be  done  by  defining  Z  (T )  ■  Z(I-T) 
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in  (1).  Z  (T )  Is  also  submodular,  and  results  from  Section  3  would  remain 
valid  for  z' .  The  reverse  algorithm  might  he  better  suited  to  handle  problems 
where  an  optimal  number  of  warehouses  tends  to  be  close  to  the  total  number 
of  potential  warehouse  sites. 

6.  Example . 

Figure  4  shows  the  cost  tableau  for  a  2x4*5  multistage  location 
problem.  Fixed  charges  for  warehouses  1-5  are:  150,  217.  200,  264,  and  140. 

For  an  explanation  of  how  the  costs  in  the  tableau  are  calculated,  see  Figure  1 
and  Section  4.  The  diagram  in  Figure  6  is  a  5 -dimensional  hypercube  whose  nodes 
represent  all  of  the  possible  open  warehouse  combinations.  The  search  tree  will 
be  a  subgraph  of  this  hypercube.  The  number  above  each  node  is  the  value  of 
an  optimal  solution  to  P  when  only  those  warehouses  marked  in  circle  of  the 
node  are  open.  The  number  in  the  parenthesis  above  a  node  is  the  total  cost 
lower  bound  for  that  node  obtained  by  calculating  (20). 


We  illustrate  the  steps  of  the  algorithm  applied  to  this  example. 
Step  Calculations 


(l) 

Set  C-l,  Q-0.  ZU  •  Z(0)  -  2107. 

Applying  Rule 

1  we 

get 

iZ^I  -  { 1 } )  •  2303  -  2201  -  102; 

iZ2(l  -  (2})  • 

217  ; 

iZ3(I  -  (3])  -  316; 

iZ4(l  -  (4 } )  -  162;  iZ5(I  -  {5}) 

•  140.  1  ■  Q 

-  0. 

For  each 

1  <  ’.1,2, 3, 4, 5}  apply  test  (20).  None  of  the  lower  bounds  exceeds  the 
current  upper  bound.  (See  the  numbers  in  parenthesis  in  Figure  6.)  Next 
calculate  iZjd)  -  1880  -  2107  •  -227;  iZ2(0)  -  -94;  ^(0)  - 

AZ4(9)  -  -136;  AZ5(9)  -  -105.  Let  List  (l)  -  (2, 5, 3, 4}.  Let 

z"  •  1880;  K  •  2 ;  go  to  Step  (2). 


-123; 

Q  •  ( 1 }  i 
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Stcp  Calculation* 

(2)  Z(I)  •  2334;  If*  1002;  the  fixed  charge  ordering  is  f  -  140, 

1 « I  1  5 

t,  •  150,  f,  ■  200,  f,  »  217,  f,  ■  264.  The  DV  ordering  is  DV_  »  0, 

1  4  2  j  2 

DVr  -  0.  DV  -  38,  DV,  •  42,  DV,  -  48.  Z°  -  1880.  Setting  k*  -  2 
?  i  3 

w*  get  t  rjs  (19)  1332  +  290  +  38  <.  1880  <  1332  ♦  490  ♦  0. 

Since  3.  •  1  2,  let  (.  •  2 ;  let  Q  •  {l};  go  to  Step  (3). 

<)/  For  each  1  <  (2, 3, 4, 3)  apply  test  (20)  setting  T  ■  (l)  U  {i'j.  (See 

Figure  6.)  For  each  1  <  {2, 5,4,3}  calculate  iZj({l}).  AZj({l})  ■  -118 

2.2 .  < { 1  i )  •  32;  4Z .  ( { 1 } )  -  22;  i2_ ( x  1  ^ )  --18.  Ut  Zu  •  1762. 

4  5  2 

List  (2)  -  {2}.  Go  to  Step  (4). 

(2)  Setting  k*  •  2  we  get  from  (19)  1332  a  290  +  38  <  1762  <  1332  +490+0 

Since  K  •  2  >2  go  to  Step  (4). 


(4) 

Let  K  ■  1 .  Remove 

(1) 

from  Q 

(now 

Q  -  0). 

Ut  k  •  4 .  Ut 

List  (1)  -  {2,5,3}. 

Ut 

Q  •  l*}. 

Ut 

-  2; 

go  to  Step  (3). 

(3)  For  each  1  «  {2,5,3}  apply  test  (2)  setting  T  “  {4}  U  {i}.  All 
solutions  are  fathomed  (see  --jfiibers  in  parenthesis  ^.n  Figure  6). 

List  (2)  •  J.  Go  to  Step  (4). 

(4)  Let  i  •  1.  Remove  {4}  from  Q  (now  Q  *  0).  Let  k  ■  3.  Let 

List  (1)  •  {2,5}.  Let  Q  -  {3).  Ut  t  •  2.  Go  to  Step  (3). 

(3)  For  each  1  <  {2,5}  apply  test  (20)  setting  T  •  (3}  U  {i}.  All 
solutions  are  fathomed  (see  Figure  6).  List  (2)  *  0.  Go  to  Step  (4). 

(4)  Ut  L  •  1 .  Remove  (3}  from  Q.  (now  Q  «  0).  Ut  k  ■  5.  Ut  List  (1)  • 

{2},  Ut  Q  •  {5}.  Ut  4*2;  Go  to  Step  (3). 
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Step  Calculations 

(3)  For  1  «  {21  apply  test  (20)  setting  T  ■  (51  U  {i}.  The  solution 
is  fathomed  (see  Figure  6).  List  (2)  »  0.  Co  to  Step  (It). 

(!t)  Let  t  ■  1.  Co  to  Step  (5). 

(5)  Stop,  z'  “  1 7 1>2  is  optimal  for  warehouses  (1.31.  An  optimal  solution 

is  shown  in  Figure  3. 

An  optimal  solution  to  the  example  is  given  In  Figure  3  when  warehouses 

1  and  3  are  opened.  The  greedy  solution  is  optimal  and  is  obtained  at  the  circled 

vertex  {1.3}  in  Figure  6.  The  upper  bound  on  an  optimal  number  of  warehouses, 

* 

k  •  2,  was  obtained  in  the  first  application  of  Step  (2)  and  thus  none  of 
level  three  warehouse  combinations  were  examined.  In  total,  16  of  the  vertices 
in  Figure  6  were  generated.  Six  of  the  vertices  on  level  2  were  fathomed  using 
(20).  The  total  number  of  transportation  pivots  required  to  solve  the  sample 
problem  was  85 . 

7 .  Computat Iona  1  Results. 

The  cost  operator  algorithm  of  Section  5  was  coded  in  FORTRAN  IV  and 
the  runs  were  made  on  a  DEC  20  time  sharing  system.  Many  of  the  problems  were 
run  at  different  times  during  the  day  and  thus  the  execution  times  may  vary  up 
to  ten  or  fifteen  percent  depending  upon  the  computing  load  of  the  machine. 

The  maximum  time  alloted  for  solving  any  problem  was  500  seconds. 

All  of  the  problems  are  derived  from  the  Kuehn  and  Hamburger  data  which 
was  originally  presented  In  [14].  Problem  sets  I  through  VII  are  taken  from 
Aklnc  and  Khumawala  [1],  and  VIII  and  IX  are  taken  from  Ellweln  and  Cray  [8]. 
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These  are  single  stage  problems  of  the  type  PI.  Problem  sets  X  through  XV 
are  multistage  problems  which  are  solved  here  for  the  first  time. 

Since  it  is  necessary  to  read  several  articles  to  determine  how  the 
problems  I  through  VII  were  originally  created,  we  will  describe  them  here. 

The  Kuehn  and  Hamburger  data  represents  a  multistage  system  with  3  factories, 

2<*  potential  warehouse  sites,  and  50  demand  centers  located  in  the  continental 
I'nlted  States.  Problems  I  through  IV  were  formed  by  considering  only  the  first 
15  potential  warehouse  sites  (l.e.,  Atlanta,  Boston,  through  New  Orleans  as  in 
l  la]).  To  obtain  the  shipping  costs  (c^)  from  these  sites  to  the  50  demand 
centers  multiply  each  distance  by  $. 025/mile  unit  cost,  which  is  the  bulk 
shipping  rate.  The  sixteenth  warehouse  is  actually  the  factory  at  Indianapolis. 
V.  used  the  distance  from  Indianapolis  to  the  50  demand  centers  multiplied  by 
S. 0125/mile  unit  as  its  shipping  cost.  The  factory  at  Indianapolis  has  the 
same  capacity  as  the  other  warehouses  and  it  nas  a  zero  fixed  charge.  The 
problems  V  through  VII  are  set  up  in  a  similar  fashion  except  that  all  2* 
potential  warehouse  sites  are  used.  Problems  X  through  XV  are  multistage 
problems  which  are  also  derived  from  the  Kuehn  and  Hamburger  data.  All  three 
factory  sites  were  used.  The  factories  have  capacities  of  30,000  or  35,000 
units.  To  derive  the  shipping  costs,  we  multiplied  the  distances  from  each 
factory  to  each  warehouse  by  $.0125/mile  unit  cost,  and  from  each  factory  to 
each  demand  center  by  .0375/mlle  unit  cost.  To  get  the  shipping  costs  from  the 
warehouses  to  the  demand  centers,  we  multiplied  the  distance  by  $.0250/mlle. 

Table  1  contains  a  description  of  all  of  the  test  problems.  Table  2 
contains  the  computational  results  for  the  single  stage  problems  and  Table  3  for 
the  multistage  problems.  The  percent  error  formula  used  to  evaluate  the  greedy 


solution  was . 
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100. 


G  * 

where  2  Is  the  greedy  value  and  Z  is  the  optimal  value.  Also  included 
in  Table  2  are  the  execution  times  tor  the  same  problems  solved  in  |1]  and  [17], 
For  test  purposes  the  problem  sets  1  through  XV  are  very  interesting  be¬ 
cause  they  contain  problems  with  a  wide  range  of  difficulty. 

The  results  for  the  greedy  solution  were  startling.  The  largest  percent 
error  Incurred  was  3.7  (see  Table  2,  IX),  however  in  general  the  percent  error 
was  less  than  .3  percent  and  the  greedy  solution  was  an  optimal  solution  in  32 
out  of  51  problems. 

Another  interesting  statistic  is  the  number  of  vertices  required  to  find 

the  optimal  solution  as  compared  to  the  total  number  of  vertices  searched.  An 

optimal  solution  is  usually  located  very  early  in  the  search  process  and  most 

of  the  effort  is  typically  spent  verifying  its  opt  lea  1 1 ty .  This  experience  is 

similar  to  that  found  in  solving  other  kinds  of  Integer  programming  problems. 

* 

The  set  I  .  which  represents  the  warehouses  fixed  open  in  the  Initializa¬ 
tion  phase  of  the  cost  operator  algorithm,  often  accounts  for  as  much  as  85^ 
of  the  total  number  of  open  warehouses  in  an  optimal  solution.  The  uncapacitated 
problems  (sets  IV.  VII,  XII,  XV)  were  very  easy  to  solve  as  expected.  As  far 
as  problem  difficulty  is  concerned,  those  problems  for  which  an  optimal  number 
of  warehouses  is  approximately  one  half  of  the  total  number  of  potential  ware¬ 
house  sites  are  usually  the  most  difficult.  This  is  probably  due  to  the  fact 
that  the  number  of  possible  warehouse  comb inat ions ,  when  t  is 

(/*■!)  ,,  when  L  is  odd,  is  the  maximum  of  the  binomial  coefficients  ^ 


even  or 


a 


The  C PC  times  quoted  in  Tables  2  were  obtained  by  three  different  sets  of 
authors  and  are  difficult  to  compare  due  to  the  differences  in  computers  and  in 
programming  efficiency.  As  far  as  machine  speeds  is  concerned,  the  IBM  370/168 
is  the  fastest,  followed  by  the  IBM  370/ 1  tj 5  and  then  by  the  DEC-20.  However, 
actual  speed  ratio  factors  for  pairs  of  the  computing  machines  are  virtually 
impossible  to  find.  It  seems  fair  to  state  that  the  performance  of  the  three 
codes  shown  in  the  table  are  not  significantly  different;  we  may  say  they 
represent  the  state  of  the  art  of  computational  results  on  these  problems. 

Table  3  exhibits  computat lonal  results  on  multistage  problems  having 
3  factories,  2-  warehouses,  and  50  demand  centers.  As  can  be  noted,  the  computa¬ 
tion  times  are  relatively  small  and  exhibit  a  relatively  small  variance  for 
integer  programming  problems.  We  again  found  the  performance  of  the  greedy 
solution  to  be  even  better  than  for  the  single  stage  problems.  This  is  perhaps 
due  to  the  fact  that  the  factories  have  a  large  total  capacity,  and  they  can 
ship  to  customers  directly  if  many  of  the  warehouses  arc  closed  down. 

8  .  Cone lu> ions . 

We  have  presented  a  cost  operator  algorithm  for  solving  multistage  loca- 

I 

t lon-a l loca t ion  problems  which  does  not  employ  problem  relaxations  as  do  the 
other  currently  best  approaches  [l.  17).  Computat iona 1  results  indicate  that 
this  method  is  competitive  with  the  others.  The  greedy  solutions  obtained  by 
the  method  are  usually  extremely  close  to  or  are  optimal.  Also,  because  the 
method  computes  many  near  optimal  solutions  as  It  solves  the  problem,  these  near 
optimal  solutions  can  be  saved  and  printed  out  for  use  by  a  manager  if  he  desires. 
Computat lonal  results  on  the  solution  of  multistage  location  problems  are  presented 
here,  but  we  have  been  unable  to  find  other  published  results  on  such  problems  for 
comparison.  The  performance  of  our  method  on  these  problems,  is  encouraging. 

We  wish  to  thank  Professor  I'mlt  Aklnc  for  supplying  the  Kuehn-Hamburger 
data  used  In  this  paper.  { 
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Table 

1 

Te*t  Problems 

Problem  Set 

*  Problems 

(q*t‘r ) 

A1 

S, 

l 

fi 

I 

4 

0« 16*50 

58268 

- 

5000 

7,500/12,500/ 

17,500/25.000/ 

II 

1 

0* 16*50 

58268 

- 

10000 

17,500 

III 

4 

0* 16' 50 

58268 

- 

15000 

7.500/12.500/ 

17.500/25,000/ 

IV 

4 

0* lb* 50 

58268 

- 

58628 

7,500/12.500/ 

17,500/25.000/ 

V 

4 

0*25' 50 

58268 

- 

5000 

7,500/12.500/ 

17.500/25,000/ 

VI 

4 

0*25*50 

58268 

- 

15000 

7.500/12,500/ 

17,500/25.000/ 

VII 

4 

0*25*50 

58268 

- 

58628 

7,500/12,500/ 

17.500/25,000/ 

VIII* 

i 

0*15*45 

3  2440 

- 

1000 

5000J 

15,000 

40,000,’ 

* 

IX 

i 

0*15*45 

37440 

- 

1500 

7500/ 

22,500 

60,000, 

X 

4 

3*24*50 

58268 

35000 

5000 

7,500/12,500/ 

17.500/25.000/ 

XI 

u 

3*24*50 

58268 

35000 

15000 

7 , 500/ 12 , 500/ 
17,500/25,000/ 

XII 

4 

3*24*50 

58268 

35000 

58628 

7,500/12.500/ 

17,500/25,000/ 

XIII 

4 

3*24*50 

58268 

30000 

5000 

7,500/12,500/ 

17,500/25,000/ 

XIV 

4 

3*24*50 

58268 

30000 

15000 

7.500/12,500/ 

17.500/25.000/ 

XV 

4 

3*24*50 

58268 

30000 

58268 

7,500/12,500/ 

17,500/25,000/ 

* 

Number*  In  the  bracket*  correspond  to  range*. 

4 


Computation  terminated 
Not  available 


Table  3 

Computet  tonal  Remits 


\  ♦  vertices  Total 

Problem  error  to  # 

Set  Greedy  Opt imal  Vertices 


ill 


Optimal 

# 

Warehouses 


Total 

# 

Pivots 


CPU 

Time 

(Dec-20) 


12345 
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